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Systematics-insensitive periodic signal search with K2 
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ABSTRACT 


From pulsating stars to transiting exoplanets, the search for periodic signals in 
K2 data, Kepler’s 2-wheeled extension, is relevant to a long list of scientihc goals. 
Systematics affecting K2 light curves due to the decreased spacecraft pointing 
precision inhibit the easy extraction of periodic signals from the data. We here 
develop a method for producing periodograms of K2 light curves that are insen¬ 
sitive to pointing-induced systematics; the Systematics-insensitive Periodogram 
(SIP). Traditional sine-htting periodograms use a generative model to hnd the 
frequency of a sinusoid that best describes the data. We extend this principle 
by including systematic trends, based on a set of ‘Eigen light curves’, following 


Foreman-Mackey et ah (2015), in our generative model as well as a sum of sine 


and cosine functions over a grid of frequencies. Using this method we are able to 
produce periodograms with vastly reduced systematic features. The quality of the 
resulting periodograms are such that we can recover acoustic oscillations in giant 
stars and measure stellar rotation periods without the need for any detrending. 
The algorithm is also applicable to the detection of other periodic phenomena 
such as variable stars, eclipsing binaries and short-period exoplanet candidates. 
The SIP code is available at https://github.coni/RuthAngus/SIPK2, 


Subject headings: asteroseismology — methods: statistical — methods: data 
analysis — stars: rotation — techniques: photometric 


1. Introduction 

The excellent precision achieved by the original Kepler mission relied on extremely 
precise pointing, for which three reaction wheels were required. After the failure of one 
of these wheels, the Kepler team devised a new pointing scheme in which the spacecraft 
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is stabilized by the Solar wind for ecliptic plane viewing zones (Howell et al. 2014). In 


this conhguration the spacecraft is able to maintain an unstable equilibrium, with the two 
functioning reaction wheels controlling pitch and yaw whilst the spacecraft slowly rolls about 
the boresight. The spacecraft hres its thrusters once every ~ 6 hours (Vanderburg and 


Johnson 2014, hereafter VJ14) to correct for this slow drift and, as stars move across pixels 
with different sensitivities, their flux varies. The extraction of high-precision photometry 
from K2 target pixel hies, despite the reduced pointing precision, is a requirement for many 
helds of research and several methods for the extraction and detrending of K2 light curves 
have already been developed. For example, VJ14 and Crossheld et al. (2015) use simple 


aperture photometry and correct the light curve of each star individually and [Aigrain et al. 


(2015) use a Gaussian process to model the non-linear dependence of stellar hux on the roll 


angle of the telescope. 

Whilst these methods successfully remove most systematic trends and produce light 
curves suitable for exoplanet search and some stellar variability studies, residual systematics 
can still affect the light curves on timescales relevent to asteroseismology and stellar rota¬ 
tion. In particular, the ~ 6 hour thruster bring signal may still appear with high power in 
the periodograms of these detrended light curves (see hgurej^. A detrending method for 
K2 light curves, specihcally intended for the asteroseismic analysis of giant stars has been 


developed by Lund et al. (2015), in which the systematics due to roll are corrected, again 


on a star-by-star basis and any remaining periodic signals at 47 /xHz (6 hour period) or 
its harmonics are removed by prewhitening. The method developed here, the Systematics- 
Insenstive Periodogram (SIP) produces periodograms of K2 light curves without the need 
for detrending or prewhitening. 


1.1. Asteroseismology 


As well as providing data that have lead to the discovery thousands of exoplanets, the 
original Kepler mission revolutionized many helds of stellar astronomy, particularly aster¬ 
oseismology. Fundamental stellar parameters—in some cases, extremely precise ones—can 
be calculated for Kepler asteroseismic stars from the power spectra of their light curves. 
Although Sun-like stars oscillate at high frequencies and require short-cadence observations, 
pulsations of giant stars he below the Nyquist frequency set by the 28.5 minute sampling 
rate of long cadence Kepler data: 283 /iHz. Asteroseismic analysis of data from the original 
Kepler mission is traditionally conducted upon detrended light curves. For short cadence 


Kepler data, this detrending method is described in Garcia et al. (2011). Due to the precise 


pointing of the original Kepler mission, systematics present in these light curves, caused by 




















3 


EPIC 201183188 



.Lj, „.-u..... 

50 100 150 200 250 


u (/iHz) 


Fig. 1.— LS periodograms of the raw (Top) and VJ14-detrended (middle) K2 light curves 
of EPIC 201183188. The bottom panel shows the SIP for this target. Peaks at ~ 47 /iHz 
and its harmonics produced by the regular spacecraft thruster hres are still present in the 
LS periodogram of the detrended data, but do not appear in the SIP. 
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temperature fluctuations and minor pointing shifts, are relatively low amplitude. 

However, this is not the case for K2 light curves: the precision over a 6 hour timescale 
is estimated to be 4 times worse in K2 data (Howell et ah [2014 ), therefore new approaches 
to the treatment of systematics are necessary. Figure [T] demonstrates the need for careful 
systematics treatment of K2 photometry for asteroseismology. The top panel shows a Lomb- 
Scargl^ (LS) periodogram of the raw, simple aperture photometry]^ of EPIC 201183188, a 
pulsating giant star. The large peaks at ~ 47 /iHz and its harmonics are caused by the 
regular thruster hres of the spacecraft. The bottom panel shows the LS periodogram of 
this light curve, after it has been detrended using the method of VJ14. The large peaks 
are still present in the detrended light curve. While this remaining noise source does not 
interfere with the detection of high-signal-to-noise transit events for periods greater than ~1 


day (Vanderburg et ah 2015), it does hamper the detection of smaller signals, particularly 


on time scales comparable to that of thruster hres. These peaks lie in an important region 
of parameter space for giant star asteroseismology and could affect the stellar parameters 
measured for thousands of giants if not dealt with appropriately. 


1.2. Stellar rotation 


Stellar rotation studies have hugely benehtted from the era of high-precision space pho¬ 
tometry. Active regions on the surface of rotating stars produce periodic variations in hux 
and stellar rotation periods can therefore be measured from Kepler light curves. Stellar 
rotation is a held of active interest as the rotation period of star can be used to infer its age 
via gyro chronology ([Skumanich 1972; Barnes 2007 Epstein and Pinsonneault 2014; Angus 


et al. 2015), is thought to be tied to the stellar magnetic dynamo, and could even reveal 


dynamical interations with companion stars or planets (e.g. Beky et al. 2014 Poppenhaeger 


and Wolk 2014). Current methods for measuring rotation periods from Kepler light curves 


include periodogram (e.g. Reinhold and Reiners 2013), Autocorrelation Function (ACF) 


( McQuillan et al.]|2013a ) and wavelet (e.g. Garcia et al.||2014 ) analysis, or some combination 
thereof. Stellar variability is not typically sinusoidal, therefore sine-htting periodograms are 
not perfectly suited to measuring rotation periods ( McQuillan et al.||2013a ). For this reason, 
the ACF method is often favoured over the periodogram method. However, because auto¬ 
correlation is performed directly on detrended light curves, and cannot be written down as a 


^All LS periodograms produced in this project were produced using the gatspy Python module: https: 
//github.com/astroML/gatspy/tree/master/gatspy/periodic 

^The method used to extract this photometry is described in (2 
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generative model, it is not possible use autocorrelation techniques on untreated K2 data. A 
quasi-periodic Gaussian process is a much better effective model for stellar variability than 
a sinusoid, however we choose to focus on the more generally applicable (and computation¬ 
ally tractable) sine-wave periodogram, leaving the Gaussian process periodogram for future 
consideration. 

In this article we focus on the examples of asteroseismology and stellar rotation, however 
many other fields of astronomy utilize periodic information in K2 light curves. These include 
studies of eclipsing binaries, variable stars, exoplanets, white dwarfs and even AGN. The 
development of tools for extracting periodic information from K2 data is essential if it is to 
be as revolutionary in time-domain astronomy as the original Kepler mission was. 

In ^ we outline the method behind the SIP. In ^ we apply the SIP to real K2 light 
curves, using some giant asteroseismic pulsators and rotating stars as test cases and then 
provide the results of some simple tests which show exactly how ‘insensitive’ the SIP is to 
systematic features. Finally, we demonstrate the SIP’s usefulness regarding other periodically 
varying objects in this section, before presenting our conclusions in ^ 


2. Method 


The method implemented in this article is an extention of the planet-search algorithm 


developed by Foreman-Mackey et al. (2015) (hereafter FM15). All targets observed by Kepler 
move on the GGD in the same way, therefore the systematics affecting each individual star’s 
light curve have shared properties. The FM15 method uses this fact by decomposing the 
light curves into a set of ‘Eigen Light Gurves’ (ELGs) using Principle Gomponent Analysis 
(PGA), which can be used to model any individual star’s light curve with very little loss 
of information. This process is similar to the method used to produce PDG-MAP data for 
the original Kepler mission (Stumpe et al. 2012; Smith et al. 2012). The resulting ELGs 


from campaign 1 can be used to model any campaign 1 K2 light curve, (campaign 0 ELGs 
for campaign 0, etc) and specifically, can model the data in combination with an arbitrary 
physical model. 

In order to construct sets of ELGs for campaigns 0 and 1, FM15 downloaded the target 
pixel files for all stars in these two fields. The position of each star was predicted using the 
World Goordinate System (WGS) and 10 circular apertures placed around the star with radii 
varying from 1 to 5 pixels in steps of 0.5 pixels. Following the procedure of VJ14, the aperture 


producing the light curve with the lowest GDPP within a 6 hour window (Ghristiansen et al. 
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2012) was 


ELCs. 


selectee^ PCA was then performed on the full set of targets in 


order to produce 


FM15 used 150 of these ELCs, plus a transit model, in order to search for exoplanet 
candidates without the need for a ‘detrending’ step. The likelihood of the data, conditioned 
on the ELC-plus-transit model was calculated over a hne grid of periods and transit depths, 
resulting in the detection of 36 new exoplanet candidates. We use a very similar technique 
to hud periodic signals in K2 data. The primary difference is that we use a sinusoid rather 
than a transit model. This model is linear, therefore the likelihood function conditioned on a 
specihe frequency can be calculated and the systematics model marginalized over analytically. 

Following the notation in FM15, our model for the fcth star can be written 

fk = Awk + noise, (1) 

where fk is the vector of N flux values, 

fk = fk,2, fk,3, fk,N)'^ (2) 


at times 


tk — 

(^1, h, ts, .. 



(3) 

A is the design matrix; 







/ Xi^i 

X2,l ■ ■ ■ 

3^150,1 1 

sin(27ri2fi) 

cos(27rz/ti) \ 


A = 

Xl,2 

X2,2 ■ ■ ■ 

3^150,2 1 

sin(27ri2f2) 

cos(27rz/f2 

(4) 


\ Xi^N 

X2,N ■ • • 

2^150,W 1 

sin(27rz/tAr) 

COs(27ri/fAr) / 



where the XijS are the ELC^ with i denoting the ELC number and j the time index. The 
design matrix contains the basis functions of the linear model. The basis functions for the 
systematic features in the light curves are the ELC values at each time index, the sine and 
cosine terms are the basis functions of the sinusoidal signal of interest, and the ‘I’s describe 
a linear offset. Any K2 light curve can be reproduced as a linear combination of these basis 
functions. We are interested in the last two elements of the weight vector: the coefficients 
of the sinusoidal signal. The maximum likelihood solution for the weight vector, w is 

Wk* ^ (A^A)-iA^fk. (5) 










7 


Under this linear model with Gaussian uncertainties, the marginalized likelihood for 
the periodic amplitude is a two-dimensional Gaussian with mean given by the last two 
elements (a) of w* and covariance given by the bottom right two-by-two block (Sa) of 
where the C matrix contains observational uncertainties on the diagonal. 
These uncertainties are estimated as 1.48x the Median Absolute Deviation (MAD), following 


Aigrain et ah (2015). Therefore, the signal-to-noise ratio, S/N of the amplitude measurement 
is \/a^Sa~^a. The (S'/A)^ is calculated over a grid of frequencies to produce an SIP. The 
(S/N) operation takes into account the goodness of £t, i.e. if the amplitude of the sinusoid 
at a given frequency is not well constrained, it is penalized. 


3. Application to real light curves 


An example LS periodogram of the raw K2 photometry for giant star, EPIG 201183188 
is shown in figure Peaks appearing at 47 /xHz and its harmonics are produced by the 
regular ~ 6 hour thruster fires that repoint the spacecraft. These peaks are also present in 
periodograms of the VJ14 detrended light curves. The presence of systematic signals at these 
timescales are problematic for asteroseismic analysis since they lie in a region of frequency 
space that is often populated by giant asteroseismic modes. It is possible to remove these 
signals by ‘prewhitening’ the data, i.e. subtracting a sinusoid of that frequency from the data, 
however this process will artificially supress all signals, both systematic and astrophysical, 
at that frequency. The SIP method eliminates the necessity for any such procedure. The 
bottom panel of figure shows the SIP for the same star, demonstrating the ability of the 
SIP method to produce periodograms that are free from thruster firing signals. 


In order to search for high signal-to-noise asteroseismic modes in the giant star candi¬ 


dates of GO1059, we searched for a power excess in the SIPs using the method of Huber 


et al.| (2009): autocorrelation functions were calculated for sections of the SIP in order to 
search for regions of increased correlation and locate the frequency of maximum power. The 
increased correlation arises from the even frequency spacing of acoustic modes, and the fre¬ 
quency of maximum correlation at the location of the power excess corresponds to the large 
frequecy separation, Au. Figures [2(a)| to 2(e) show example power spectra of 6 targets for 
which we detect pulsations using this method. 


The top panel of figure shows the VJ14-detrended light curve of an active, rotating 
star, EPIG 201133037, with a linear trend subtracted off. The brightness fluctuations clearly 
visible in the light curve of this target are produced by cool active regions on the stellar 
surface, which reduce the stellar flux periodically. The rotation period of this star is therefore 
around 20 days. The middle and bottom panels show an AGE and LS periodogram of the 
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Fig. 2.— SIPs of 6 long cadence K2 giants with asteroseismic oscillations. These were 
selected from the guest observing program, GO1059 and identihed using the method of 


Huber et ah (2009) 
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detrended light curve. The top panel of hgure|^ shows the raw light curve of the same target 
in grey, with the conditioned light curve in black. This conditioned light curve was produced 
by removing the best htting systematic trends, described by a certain combination of the 
ELCs, at the best htting period of the sinusoid. The bottom panel shows the SIP. 

Each of these three methods measures a rotation period of around 20 days for this target. 
This example demonstates the ability of the SIP to recover rotation periods that agree with 
those measured from detrended light curves by autocorrelation. We also include an example 
that demonstrates the ability of the SIP to outperform a periodogram of detrended data. 
Figure shows the light curve, ACE and LS periodogram of another rotating star, EPIC 
201142043 and hgure[^ shows its SIP. This star shows lower amplitude variability than the 
previous example and the careful treatment of systematics is much more important. Whereas 
the ACE method is able to measure a rotation of ~ 3 days for this star, the LS periodogram 
of the detrended light curve incorrectly measures a period of 59 days. Although there is a 
small peak at the rotation period of the star, it is not the dominant periodic signal. The SIP 
method is, by dehnition, insensitive to these long-term systematics and is able to measure a 
period of ~ 2 days. This example further demonstrates the fact that long-term systematic 
trends caused by slow pointing variations are often not removed by conventional detrending 
methods. The 59 day signal is almost certainly a systematic trend and not an astrophysical 
signal because it does not appear in the SIP. It is well described by the ELCs and must 
therefore be common to many stars. 


We have shown that the SIP method is able to measure stellar rotation periods and 
does better than producing periodograms from detrended data. However, it has been shown 
that the ACE method often performs better than periodogram methods in general for mea¬ 
suring stellar rotation periods (McQuillan et al. 2013a||b Mazeh et al. 2015). For stars with 
relatively high-amplitude variability, for which perfect removal of systematic trends is less 
important, performing the ACE method on detrended data is likely to produce similar re¬ 
sults to the SIP method. The SIP method is ideally suited to low-amplitude cases, where 
systematic trends could drastically influence rotation period measurements. Whilst the SIP 
method may outperform ACE in the low-amplitude cases, any ‘marginal’ rotation period 
measurements calculated using either method should be treated with caution unless a rep¬ 
resentative uncertainty is provided. In general neither ACE nor periodogram methods are 
equipped to provide such uncertainties. In practise, we recommend using both the SIP and 
ACE methods, in combination with a by-eye check, to measure rotation periods for K2 stars. 
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Period (days) 

Fig. 3.— Top: Light curve of EPIC 201133037, detrended using the method of VJ14. Middle: 
Autocorrelation function of the detrended light curve. The autocorrelation function method 
measures a rotation period of 21 days for this star. Bottom: The LS periodogram of the 
detrended light curve. The highest peak in the periodogram is located at 25 days. 






11 


cu 

_> 


40000 

39800 

39600 

39400 

39200 

39000 

38800 

38600 

38400 


EPIC 201133037 

1980 1990 2000 2010 2020 2030 2040 2050 



Period (days) 


Fig. 4.— Top-. The raw light curve of EPIC 201133037 is shown in grey and the conditioned 
light curve is shown in black. The conditioned light curve is produced by removing the 
trends that best describe the data, at the best htting frequency. Bottom-. An SIP of the raw 
light curve, produced by modelling the data using the top 150 ELCs plus a sine and cosine 
function at a range of frequencies. The highest peak in the SIP is located at 25 days. 
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Fig. 5.— Top: Light curve of EPIC 201142043, detrended using the method of VJ14. Middle: 
Autocorrelation function of the detrended light curve. The autocorrelation function method 
measures a rotation period of 3 days for this star. Bottom: The LS periodogram of the 
detrended light curve. The highest peak in the periodogram is located at 59 days and is 
likely to be a systematic trend produced by spacecraft pointing variations. 
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Fig. 6.— Top-. The raw light curve of EPIC 201142043 is shown in grey and the conditioned 
light curve is shown in black. The conditioned light curve is produced by removing the 
trends that best describe the data, at the best htting frequency. Bottom: An SIP, produced 
by modelling the data as a linear combination of the top 150 ELCs plus a sine and cosine 
function at a range of frequencies, measuring a rotation period of 2 days. The SIP is, by 
dehnition, insensitive to the long-timescale systematics that dominate the LS periodogram 
of the detrended data, shown in hgure 
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3.1. Tests and discussion 


In order to demonstrate the consistent ability of the SIP method to remove the signal at 
47/rHz, corresponding to the periodic ~6 hour thruster brings, we computed SIPs for 4923 
targets from the GO1049 proposal: “Galactic Archaeology on a grand scale” (PI: Stello, D.). 
For each target, an SIP of its raw photometry and a LS periodogram of its VJ14 light curve 
was calculated for frequencies between 40 and 54 pHz. Both the height and frequency of 
the highest peak in the SIP and the highest peak in the LS periodogram were recorded. A 
histogram of the frequencies of the highest peaks in the SIPs of all 4923 targets is shown in 
the top panel of bgure The bottom panel shows the histograms of peak heights within 
the correspondingly colored ranges indicated in the top panel. This bgure shows that while 
there are a greater number of maximum peaks around 47 /iHz, the S/Ns of these peaks are 
comparable to those found just above and just below this frequency. Figure shows the 
equivalent results for the VJ14 light curves. There is a signibcant number of large peaks 
at ~47 /xHz in the LS periodograms of the detrended light curves; the highest peak in the 
LS periodograms was almost always located at ~ 47 /iHz. Furthermore, the distribution of 
peak power within the range 46.5-48 /iHz is skewed towards higher powers, i.e. a substantial 
fraction of the peaks at ~ 47 /iHz have a large power. The SIP method is able to consistently 
remove the 47 /iHz signal which is present in almost every VJ14 light curve. 


Figures and show the conditioned light curves and SIPs of an RR Lyrae star, 
selected from Guest observer program GO1018; “Long Gadence RR Lyrae targets” (PL 
Plachy, E.) and an Eclipsing Binary (EB), selected from the Armstrong et ah (2015) K2 EB 
and variable star catalogue respectively. Figure shows the conditioned light curve and 


SIP of the short-period, disintegrating planet candidate discovered by Sanchis-Ojeda et ah 


(2015). This planet has a period of around 9 hours, short enough to be detectable with a 
periodogram, as was demonstrated for a number of ultra-short period Kepler exoplanets by 


Sanchis-Ojeda et ah (2014). The top panels of these three bgures show the K2 light curves of 


these objects, conditioned on the highest S/N sinusoidal signal in the periodograms. These 
light curves were produced by subtracting the trends that best describe the data at the 
highest S/N period found in the SIP. 


Photometric variability in dwarf stars on timescales less than 8 hours, often known 
as bicker, has been linked to surface gravity ([Bastien et al. 2013 Kipping et al. 2014). 


The metrics used to quantify photometric variability include bnding the range in intensity, 
counting the number of zero crossings and calculating the root-mean-square (RMS) of the 
light curve. Although these features are related to signal processing, they are operations 
performed on detrended light curves, not inferred from periodograms. However, it may be 
possible to derive a property of the periodogram that scales with the density or surface 
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gravity of a star, for example, the mean excess power at frequencies near those relevent to 
granulation timescales. The SIP method presented here would be useful for such a technique. 
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ln(Maximum peak height [Relative (S/N)^]) 


Fig. 7.— Top: Histogram of the frequencies of the highest peaks in the SIPs of 4923 K2 
targets within the range 40-54 /iHz. Bottom: Histograms of peak heights within the cor¬ 
respondingly colored ranges indicated in the top panel. Whilst there is a larger number 
maximum peaks around 47 /iHz (the frequency corresponding to the 6 hour thruster hre) 
the amplitudes of these maximum peaks are comparable to the maximum peak heights just 
above and just below this frequency. 
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Fig. 8.— Top: Histogram of the frequencies of the highest peaks in the LS periodograms 
of the Vanderburg and Johnson (2014) light curves of 4923 K2 targets within the range 
40-54 /xHz. Bottom: Histograms of peak heights within the correspondingly colored ranges 
indicated in the top panel. The frequency of maximum peak height was ~ 47 /xHz in almost 
every periodogram. Furthermore, the distribution of maximum peak height within the range 
46.5-48 /rHz is skewed towards higher powers, i.e. a large fraction of the peaks at ~ 47 /xHz 
have a large power. 
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Fig. 9.— Top: The light curve of RR Lyrae star, EPIC 201339783, conditioned on the 
highest amplitude sinusoidal signal found in the SIP. Bottom: The systematic-insensitive 
periodogram of this light curve. This target was selected from GO1018 (PI; Plachy, E.). 
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Fig. 10.— Top: The light curve of eclipsing binary, EPIC 201473612, conditioned on the 
highest amplitude sinusoidal signal found in the SIP. Bottom: The systematic-insensitive 
periodogram of this light curve. This target was selected from the catalogue of EBs and 


variable stars of Armstrong et ah (2015) 
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Fig. 11.— Top: The light curve of exoplanet candidate, EPIC 201637175, conditioned on 
the highest amplitude sinusoidal signal found in the periodogram. Bottom: The systematic- 
insensitive periodogram of this target. 
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4. Conclusions 

We demonstrate that modelling campaign 1 K2 photometry as a linear combination of 
150 PCA components plus a sinusoid can produce periodograms that are almost completely 
free from instrumental systematic signals, without the need for detrending. We hnd that the 
47 /iHz signal, generated by the spacecraft thruster hres is not present in the vast majority 
of Systematics-Insensitive Periodograms (SIPs) for more than 4000 targets selected from the 
K2 guest observer program, GO1059, “Galactic Archaeology on a grand scale” (PI; Stello, 
D.). The SIP method is highly successful for campaign 1 targets where the large number of 
stars, observed for a baseline of 80 days ensures that most of the systematics are captured in 
the ELGs and we anticipate that it will be equally effective for the up-and-coming campaigns. 

The SIP method is capable of detecting periodicities in K2 data in the region of frequency 
space relevent to the study of asteroseismic oscillations in giant stars and for any signals with 
a timescale close to 6 hours. It is also effective at measuring stellar rotation periods and 
is an improvement upon a simple LS periodogram of detrended data. In practise, the best 
approach for measuring rotation periods in K2 data is likely to be a combination of the 
SIP method and the AGP method, where autocorrelation is performed on detrended light 
curves. The SIP code is available for public use and can be found at https://github.coni/ 
RuthAngus/SIPK2, 

It is a pleasure to thank Dan Huber (Sydney) who provided many excellent comments for 
this paper and useful asteroseismology tips. We would also like to thank David Hogg (NYU) 
and Suzanne Aigrain (Oxford) for their comments, plus Andrew Vanderburg (Harvard), Ben 
Montet (Galtech) and Stephanie Douglas (Golumbia) for their extremely helpful suggestions 
and recommendations regarding this project. J.A.J is supported by generous grants from 
the David and Lucile Packard and Alfred P. Sloan Foundations. The data presented in 
this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is 
operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. Support for MAST for non-HST data is provided by the NASA Office 
of Space Science via grant NNX09AF08G and by other grants and contracts. This paper 
includes data collected by the Kepler mission. Funding for the Kepler mission is provided 
by the NASA Science Mission directorate. 


REFERENCES 


S. Aigrain, S. T. Hodgkin, M. J. Irwin, J. R. Lewis, and S. J. Roberts. Precise time series 
photometry for the Kepler-2.0 mission. MNRAS, 447:2880-2893, March 2015. doi: 



10.1093/mnras/stu2638. 


R. Angus, S. Aigrain, D. Foreman-Mackey, and A. McQuillan. Calibrating gyro chronology 

using Kepler asteroseismic targets. MNRAS, 450:1787-1798, June 2015. doi: 10.1093/ 
mnras/stv423. 

D. J. Armstrong, J. Kirk, K. W. F. Lam, et al. K2 Variable Catalogue II; Variable Stars 
and Eclipsing Binaries in K2 Fields 1 and 0. ArXiv e-prints, February 2015. 

S. A. Barnes. Ages for Illustrative Field Stars Using Gyrochronology: Viability, Limitations, 

and Errors. ApJ, 669:1167-1189, November 2007. doi: 10.1086/519295. 

F. A. Bastien, K. G. Stassun, G. Basri, and J. Pepper. An observational correlation between 
stellar brightness variations and surface gravity. Nature, 500:427-430, August 2013. 
doi; 10.1038/naturel2419. 

B. Beky, M. J. Holman, D. M. Kipping, and R. W. Noyes. Stellar Rotation-Planetary Orbit 

Period Commensurability in the HAT-P-11 System. ApJ, 788:1, June 2014. doi: 
10.1088/0004-637X/788/1/1. 

J. L. Christiansen, J. M. Jenkins, D. A. Caldwell, et al. The Derivation, Properties, and 
Value of Kepler” s Combined Differential Photometric Precision. PASP, 124:1279- 
1287, December 2012. doi: 10.1086/668847. 

1. J. M. Crossheld, E. Petigura, J. Schlieder, et al. A nearby M star with three transiting 
super-Earths discovered by K2. ArXiv e-prints, January 2015. 

C. R. Epstein and M. H. Pinsonneault. How Good a Clock is Rotation? The Stellar Rotation- 

Mass-Age Relationship for Old Field Stars. ApJ, 780:159, January 2014. doi; 10.1088/ 
0004-637X/780/2/159. 

D. Foreman-Mackey, B. T. Montet, D. W. Hogg, et al. A systematic search for transiting 

planets in the K2 data. ArXiv e-prints, February 2015. 

R. A. Garcia, S. Hekker, D. Stello, J. Gutierrez-Soto, R. Handberg, D. Huber, G. Karoff, 
K. Uytterhoeven, T. Appourchaux, W. J. Ghaplin, Y. Elsworth, S. Mathur, J. Ballot, 
J. Ghristensen-Dalsgaard, R. L. Gilliland, G. Houdek, J. M. Jenkins, H. Kjeldsen, 
S. McGauliff, T. Metcalfe, G. K. Middour, J. Molenda-Zakowicz, M. J. P. F. G. 
Monteiro, J. G. Smith, and M. J. Thompson. Preparation of Kepler light curves for 
asteroseismic analyses. MNRAS, 414:L6-L10, June 2011. doi: 10.1111/j.1745-3933. 
2011.01042.x. 



- 23 - 

R. A. Garcia, T. Ceillier, D. Salabert, et al. Rotation and magnetism of Kepler pulsating 

solar-like stars. Towards asteroseismically calibrated age-rotation relations. A&A, 
572:A34, December 2014. doi; 10.1051/0004-6361/201423888. 

S. B. Howell, C. Sobeck, M. Haas, M. Still, T. Barclay, F. Mullally, J. Troeltzsch, S. Aigrain, 

S. T. Bryson, D. Caldwell, W. J. Chaplin, W. D. Cochran, D. Huber, G. W. Marcy, 
A. Miglio, J. R. Najita, M. Smith, J. D. Twicken, and J. J. Fortney. The K2 Mission: 
Characterization and Early Results. PASP, 126:398-408, April 2014. doi: 10.1086/ 
676406. 

D. Huber, D. Stello, T. R. Bedding, et al. Automated extraction of oscillation parameters for 
Kepler observations of solar-type stars. Communications in A stew seismology, 160: 
74, October 2009. 

D. M. Kipping, F. A. Bastien, K. G. Stassun, et al. Flicker as a Tool for Characterizing 
Planets Through Asterodensity Prohling. ApJ, 785:L32, April 2014. doi: 10.1088/ 
2041-8205/785/2/L32. 

M. N. Lund, R. Handberg, G. R. Davies, W. J. Chaplin, and C. D. Jones. K2P$"2$ $-$ A 
photometry pipeline for the K2 mission. ArXiv e-prints, April 2015. 

T. Mazeh, H. B. Perets, A. McQuillan, and E. S. Goldstein. Photometric Amplitude 

Distribution of Stellar Rotation of KOIs—Indication for Spin-Orbit Alignment of 
Cool Stars and High Obliquity for Hot Stars. ApJ, 801:3, March 2015. doi: 
10.1088/0004-637X/801/1/3. 

A. McQuillan, S. Aigrain, and T. Mazeh. Measuring the rotation period distribution of held 
M dwarfs with Kepler. MNRAS, 432:1203-1216, June 2013a. doi: 10.1093/mnras/ 
stt536. 

A. McQuillan, T. Mazeh, and S. Aigrain. Stellar Rotation Periods of the Kepler Objects of 
Interest: A Dearth of Close-in Planets around Fast Rotators. ApJ, 775:L11, September 
2013b. doi: 10.1088/2041-8205/775/1/Lll. 

K. Poppenhaeger and S. J. Wolk. Indications for an inhuence of hot Jupiters on the rotation 
and activity of their host stars. A&A, 565:L1, May 2014. doi: 10.1051/0004-6361/ 
201423454. 

T. Reinhold and A. Reiners. Fast and reliable method for measuring stellar differential 
rotation from photometric data. A&A, 557:A11, September 2013. doi: 10.1051/ 
0004-6361/201321161. 



-24- 

R. Sanchis-Ojeda, S. Rappaport, J. N. Winn, et al. A Study of the Shortest-period Planets 
Found with Kepler. ApJ, 787:47, May 2014. doi; 10.1088/0004-637X/787/1/47. 

R. Sanchis-Ojeda, S. Rappaport, E. Palle, et al. The K2-ESPRINT Project I; Discovery of 
the Disintegrating Rocky Planet with a Cometary Head and Tail EPIC 201637175b. 
ArXiv e-prints, April 2015. 

A. Skumanich. Time Scales for CA II Emission Decay, Rotational Braking, and Lithium 
Depletion. ApJ, 171:565, Eebruary 1972. doi: 10.1086/151310. 

J. C. Smith, M. C. Stumpe, J. E. Van Cleve, et al. Kepler Presearch Data Conditioning 
II - A Bayesian Approach to Systematic Error Correction. PASP, 124:1000-1014, 
September 2012. doi: 10.1086/667697. 

M. C. Stumpe, J. C. Smith, J. E. Van Cleve, et al. Kepler Presearch Data Conditioning I 
- Architecture and Algorithms for Error Correction in Kepler Light Curves. PASP, 
124:985-999, September 2012. doi: 10.1086/667698. 

A. Vanderburg and J. A. Johnson. A Technique for Extracting Highly Precise Photometry 
for the Two-Wheeled Kepler Mission. PASP, 126:948-958, December 2014. doi: 
10.1086/678764. 

A. Vanderburg, B. T. Montet, J. A. Johnson, et al. Characterizing K2 Planet Discoveries: A 
Super-Earth Transiting the Bright K Dwarf HIP 116454. ApJ, 800:59, February 2015. 
doi: 10.1088/0004-637X/800/1/59. 


This preprint was prepared with the AAS MLX macros v5.2. 



